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Abstract 



o 

"q ■ Short title: sl(2|l) Meixner oscillator 
Q '■ PACS numbers: 03.67.Hk, 02.30.Gp 

We explore a model for the one-dimensional quantum oscillator based upon the Lie superal- 
gebra st(2|l). For this purpose, a class of discrete series representations of sl(2|l) is constructed, 
each representation characterized by a real number /3 > 0. In this model, the position and mo- 
^ (— i mentum operators of the oscillator are odd elements of st(2|l) and their expressions involve an 

arbitrary parameter 7. In each representation, the spectrum of the Hamiltonian is the same as 
that of the canonical oscillator. The spectrum of the momentum operator can be continuous or 
infinite discrete, depending on the value of 7. We determine the position wavefunctions both in 
the continuous and discrete case, and discuss their properties. In the discrete case, these wave- 
£Nl " functions are given in terms of Meixner polynomials. From the embedding osp(l|2) C sl(2|l), it 

can be seen why the case 7 = 1 corresponds to the paraboson oscillator. Consequently, taking 
■ the values (^,7) = (1/2, 1) in the bI(2 1 1 ) model yields the canonical oscillator. 

in 

1 Introduction 

OCT 

o. 

\ The quantum harmonic oscillator is one of the main examples in physics, both because of its 
use in physical models and approximations and because of its mathematical beauty. Besides the 
canonical oscillator there are many non-canonical versions that received attention in the literature. 
In particular, there are many algebraic constructions to model a quantum oscillator by extending (or 
deforming) the common oscillator Lie algebra. The difficulty for such models is often to determine 



the spectra of observables and an explicit form of their eigenfunctions. Only for some models, one 
can develop such a complete theory. One of these models is the g-oscillator, a (/-deformation of the 
standard quantum oscillator |li p 29 t l37]. discussed in this context in [2] and |24| . 

During the last decade, new oscillator models were developed such that the same dynamics as 
in the classical or quantum case is satisfied, and in such a way that the operators corresponding to 
position, momentum and Hamiltonian are elements of some algebra different from the traditional 
Heisenberg (or oscillator) Lie algebra. In the one-dimensional case, there are three (essentially self- 
adjoint) operators involved: the position operator q, its corresponding momentum operator p and 
the Hamiltonian H which is the generator of time evolution. The main requirement is that these 
operators should satisfy the Hamilton-Lie equations (or the compatibility of Hamilton's equations 
with the Heisenberg equations): 

[H,q] = -ip, [H,p]=iq, (1) 



1 Permanent address: Institute of Physics, Azerbaijan National Academy of Sciences, Javid av. 33, AZ-1143 Baku, 
Azerbaijan 



1 



in units with mass and frequency both equal to 1, and h = 1. Contrary to the canonical case, the 
commutator [q,p] = i is not required. Apart from ([1]) and the self-adjointness, it is then common 
to require the following conditions [3]: 

• all operators q, p, H belong to some Lie algebra or Lie super algebra A; 

• the spectrum of H in (unitary) representations of A is equidistant. 

The model that has received most attention occurs for A = su(2) (or its enveloping algebra) [3J- 
[6]. For this case, the appropriate representations are the common su(2) representations labeled 
by an integer or half-integer j. Since these representations are finite-dimensional, one is dealing 
with "finite oscillator models", of potential use in optical image processing [6]. Recently, this 
model has been extended by introducing an additional parameter in the algebra, thus leading to 
finite paraboson oscillator models [19, 20]. In [22], the Lie algebra su(2) was extended to the Lie 
superalgebra s[(2|l), yielding another interesting finite oscillator model with appealing position 
spectra and remarkable discrete wavefunctions. 

All the examples just mentioned concern finite oscillator models, with a finite (equidistant) 
spectrum for the Hamiltonian and with a discrete spectrum for the position and momentum op- 
erator, due to the fact that the relevant (unitary) representations of the model algebra A are 
finite-dimensional. But there are also some models of quantum oscillators with continuous spectra 
of position and momentum operators available besides the canonical model. One model is based on 
the positive discrete series representations of su(l, 1) [22]. In such a representation the spectrum of 
the position operator is M. The position wavefunctions are given by normalized Meixner-Pollaczek 
polynomials. Klimyk discussed many fascinating properties of these su(l, 1) oscillators [25J. A 
^-deformation of Klimyk's su(l,l) model was investigated in [7j. The position and momentum 
operators have spectra covered by a finite interval of the real line, which depends on the value of q, 
and the wavefunctions are given in terms of g-Meixner-Pollaczek polynomials. Another extension 
of the Lie algebra su(l, 1) by means of a parity or reflection operator R was studied in [21J, together 
with a class of discrete series representations. In the corresponding model, the Hamiltonian has a 
discrete but infinite equidistant spectrum, and the position operator has spectrum E. Again, the 
model is sufficiently simple to construct the position wavefunctions explicitly, and these are given 
in terms of continuous dual Hahn polynomials |17U26j . 

In the current paper we return to the Lie superalgebra A = sl(2|l), but we shall consider a new 
class of infinite-dimensional discrete series representations labeled by a positive real number f3. In 
the model using these representations, the choice for the Hamiltonian H is unambiguous and (as 
required) it has a discrete but infinite equidistant spectrum. The choice for the position operator 
q is not completely fixed: an arbitrary real parameter 7 remains in the expression for q in terms of 
the sl(2|l) generators. The momentum operator p follows from the first equation of ([1]). Our main 
work is then devoted to determining the spectrum of q and of its formal eigenvectors from which the 
position wavefunctions ^^'^(x) follow. Our analysis shows that q has an infinite discrete spectrum 
when I7I 7^ 1, but a continuous spectrum when I7I = 1. It is quite remarkable to have both of these 
situations occurring in the same oscillator model. In the case 7 = 1, the position wavefunctions 
coincide with those of the paraboson oscillator (i.e. they are given in terms of generalized Laguerre 
polynomials), and in particular for f3 = 1/2 the model coincides with the canonical oscillator. When 
I7I 7^ 1, the position wavefunctions are given in terms of Meixner polynomials, a class of orthogonal 
polynomials with a discrete orthogonality relation. 

The structure of the paper is as follows. In the next section, we give the Lie superalgebra 
sl(2|l) and show that it possesses a class of discrete series representations. In section 3, the s[(2|l) 
oscillator model is presented. In particular, using the technique of unbounded Jacobi matrices, the 
spectrum of the position and momentum operator is determined. Section 4 is devoted to studying 
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the shape of the position wavefunctions. We examine some plots of wavefunctions, and investigate 
how the discrete wavefunctions for 7 / 1 (given in terms of Meixner polynomials) tend to the 
continuous wavefunctions for 7 = 1 (given in terms of Laguerre polynomials), both in plots and 
as a limit computation. In section 5, we determine the corresponding sl(2|l) Fourier transform, 
which is defined as the kernel relating position eigenvectors to momentum eigenvectors. Due to the 
fact that bilinear generating functions are known for the Meixner polynomials appearing here, the 
st(2|l) Fourier transform can be given in explicit form. In section 6 we briefly recall the paraboson 
oscillator model and its relation to the Lie algebra osp(l|2). The embedding osp(l|2) C sl(l|2) 
explains algebraically why 7 = 1 is a special case in our model and why it corresponds to the 
paraboson oscillator. Finally, some remarks and a further discussion of some interesting quantities 
in the model is presented in section 7. 

2 The Lie superalgebra s[(2|l) and a class of discrete series repre- 
sentations 

The Lie superalgebra st(2| 1) with even part sl(2) © is well known, and has been constructed 
in a previous paper |22| where a class of finite-dimensional representations were used. Let us just 
recall the basis here, following the choice of |13|. p. 261], consisting of four odd (or 'fermionic') basis 
elements F + , F~ ,G + ,G~ and four even (or 'bosonic') basis elements H, E + , E~ , Z , given by 

F + = e 32 , G + = e 13 , F~ = e 31 , G~ = e 23 , (2) 

1 1 

H = -(en - e 2 2), E + = e i2 , E~ = e 2 i, Z = -(en + e 22 ) + e 33 (3) 

in terms (graded) 3x3 Weyl matrices ey. The basis for the sl(2) subalgebra is {H, E + ,E~} and 
the £(t(l) = U(l) subalgebra is spanned by Z. The basic Lie superalgebra brackets can be found 
in [13 p. 261], [35] or [22|: 

{F ± ,G ± } = E ± , {F ± ,G T } = ZtH, 

{F ± ,F ± } = {G ± ,G ± } = {F±,FT} = {G ± ,G^} = 0; (4) 
[H,E ± ] = ±E ± , [E+,E-] = 2H, [Z, H] = [Z, E±] = 0; (5) 

[H,F ± ] = ± l -F ± , [Z,F ± ] = ±F ± , [£±,^ = 0, [E^,F ± ] = -F^, 

[H,G ± ] = ± l -G ± , [Z,G ± ]=-ic ± , [E ± ,G ± } = ^ [^ }G ±] = G fT. (6 ) 

The finite-dimensional irreducible representations of s[(2|l) have been studied by Scheunert et 
al [35] and Marcu [30] . but little seems to be known about infinite-dimensional representations. 
Here, we present a class of infinite-dimensional representations labeled by a positive number j3 > 0. 
We shall call them positive discrete series representations of sl(2|l), as they are closely related to 
positive discrete series representations of 5u(l, 1). 

First of all, let us fix a *-structure (or an adjoint operation) on the Lie superalgebra by 

Z+ = Z, H^ = H, ( J B ± ) t = -^ T , (F ± )t = T G T , (G ± ) t = ±F T . (7) 

It is easy to see that this ^-structure is indeed compatible with the Lie superalgebra brackets ([H)-©. 
Note that the ^-structure on sl(2) implies that we are dealing with the form su(l, 1). 

The positive discrete series representations of sl(2|l) are unitary representations labeled by /3 > 

0. The representation space is £ 2 (Z + ) equipped with an orthonormal basis |/3, n) (n = 0, 1,2, . . .), 

1. e. 

{P,m\P,n) = <W (8) 
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For the actions of the sl(2|l) basis elements on these vectors, it is handy to use the following "even" 
and "odd" functions, defined on integers n: 

£(n) = 1 if n is even and otherwise, 

0(n) = 1 if n is odd and otherwise. (9) 

Note that 0(n) = 1 — £ (n), but it is convenient to use both notations. The actions of the odd 
generators are now given by: 

F + \P,n) =8{n)^P + ^ |/3,n + l), F~\P, n) = £ (n) J| \P, n - 1), 

/ 71 -\- 1 I 71 — 1 

G+\p,n)=0(n) ] j^- \P,n + l), G~\P,n) = -0{n)^j p + — \P,n- 1). (10) 

The actions of the even generators can in principle be computed from (j4]) , and are 

Z\P,n) = -£{n)^ n) - 0(n)(^±) |/3,n), 
H\p,n) = ±(n + P) |/3, n), 



+ + |/3,n + 2) + O(„)J09 + 2±l)(!L±l) |/3,n + 2), 



/ 71 71 / 77 — 1 77 — 1 

Er|An> = -£(n)y 09 - 1 + f )(~) |/?,n- 2) - 0(n)^j{p + — )(— ) |/3,n - 2). (11) 

Proposition 1 TTie representation Hp ((3 > 0) in £ 2 (Z+) ; defined by the actions (|10[ )-( fTT1) . is an 

irreducible star representation (or unitary representation) of the Lie superalgebra sl(2|l). With 
respect to the even subalgebra sl(2) = su(l,l), it decomposes into the direct sum of two positive 
discrete series representations tip and n p+i , with Bargmann index (or label) resp. § and ^fci. 

2 2 

Proof. To show that the actions (|10p -(|li p define indeed a representation of s[(2|l) is straightfor- 
ward but tedious. Essentially, one should verify that all bracket relations (jH)-© are satisfied on the 
basis vectors \/3,n). In practice, note that (|6|) and © follow from @ and the Jacobi-identity. So 
it is sufficient to verify that (|3J) is satisfied for the action (fTUj) . and that (fTUj) implies (fTT|) using Q. 
Note that |/3, 0) is a generating vector for the representation, since 



(G+F+n/3,0) = ^m^\f3,2n), F+(G + F + ) n \p, 0) = y/n\(/3) n+1 |/3, 2n + 1), (12) 

where (/3) n is the Pochhammer symbol |1,8,36]: (P) n = /3(/3 + 1) • • • (/3 + n — 1). Irreducibility then 
follows from the actions 



(G-F-) n \P,2n) = (-l) n V^G^ I AO), (G~F~) n G~\P,2n + 1) = (-l) n+ V™K/S)n+i 1/3,0). 

(13) 

To see that the representation is a star representation for the ^-structure (J7|), it is sufficient to 
check 

(P,2n + l\F+\/3,2n) = -(/3, 2n|G"|/3, 2n + 1), 

{P,2n-l\F-\P,2n) = (/3, 2n|G + |/3, 2n - 1). (14) 
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Finally, relabel the even and odd vectors of the representation by e m = \/3,2m) and f m = 
2m + 1). Then the vectors e m (m = 0,1,2,...) are an orthonormal basis for the action of 
su(l,l) in £ 2 (Z + ), with 



He m = (m + -)e m , E + e m = ^ (0 + m)(m + 1) e m +i, E e m = -y 7 {fi + m - l)me m -i, 

so this is the representation tt§_. Note that the action of Z on this representation is — ^ times 

the identity operator. Similarly, the vectors f m (m = 0, 1, 2, . . .) are an orthonormal basis for the 
action of su(l, 1) in ^ 2 (Z+), with in particular Hf m = (m + ^^)f m , so this is the representation 
7T/3+1. In this case, the action of Z on this representation is — times the identity operator. □ 



3 An s((2|l) oscillator model and the spectrum of a position oper- 
ator 



In order to use the discrete series representations of st(2|l) for an oscillator model, it is natural to 
take the Hamiltonian H as 

H = 2H + 1 - 0. (15) 



2H + --8. 
2 

This operator is diagonal, self-adjoint, and has the equidistant spectrum: n + | (n = 0, 1, 2, . . .). 
Following the arguments of [22], we should take for the position operator q an arbitrary odd (real) 
self-adjoint element of st(2| 1) , i.e. an element of the form 



F + + 7 G+ - G~ + jF~ , 



(16) 



with 7 a real constant (an overall constant does not play a crucial role, so that is why we have 
taken the coefficient of F + equal to 1, and (|16p still represents the most general case). 
When q is fixed by (I16p . the expression of the momentum operator p follows from (pQ): 



P 



i(F+ + 7 G+ + G~ -jF-). 



(17) 



These operators (fT5|) . (fT6|) and pT|) do indeed satisfy ([T]) and the conditions described in section 1 
are satisfied, and thus we are dealing with models for the oscillator in a class of infinite dimensional 
representations of st(2|l). 

In the (ordered) basis {\0,n), n = 0,1,2,...}, the operator q is represented by an infinite 
symmetric tridiagonal matrix M q : 



( o 

Rq 



M„ 



Rq 


5*i 



Si 


Ri 



R-2 



S-2 



S 2 





(18) 



where 



Rr. 



Sn 



n 



in 



0,1,2,...). 



(19) 



For 7 > 0, such a matrix is a Jacobi matrix, and its spectral theory is related to orthogonal 
polynomials [10, 27, 28J (7 < is similar to 7 > 0: it will soon be clear that only I7I plays a role; for 
7 = the matrix (fT8j) decomposes and also that case will be easy to treat). Following the procedure 
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described in [28, s §2], one should construct polynomials p n (x) of degree n in x, with p_i(x) = 0, 
Po(x) = 1, and 

XP2n(x) = S n p 2n -l(x) + R n P2n+l (x) , 

Xp 2n +l(x) = RnP2n(x) + S n+1 p 2n +2(x), (n = 0,1,2, .. .). (20) 

Such polynomials are orthogonal for some positive weight function w(x). Then the spectrum of M q 
(or of q) is the support of this weight function. This technique works provided the (Hamburger) 
moment problem for the Jacobi matrix is determinate [101127] . This is equivalent to saying that the 
corresponding Jacobi operator is essentially self-adjoint. A sufficient condition is that |10[I27] 

x j x i x i x i 

^ n K + ^ ~S~n = ^ n 7w+^ + ^ tT^ = 00 ' 

n=0 n=l n=0 v n=l 

which is satisfied here. So the spectrum of the position operator q is just the support of the weight 
function w(x). Furthermore, for a real value x belonging to this support, the corresponding formal 
eigenvector of q is given by 



(x) = J2pn{x)\P,n). (21) 



v v 

n=0 



So the purpose is first to construct the polynomials p n (x), and then to find the corresponding 
weight function. The solution of ([20]) is given in terms of termination hyper geometric series; for 
their notation we follow that of standard books 



Proposition 2 When 7 2 7^ 1, the solution of the recurrence relations ([20]) is given by 



P 2 n(x) = (-7)-"V^^(" n '*;l 



V-7 2 



p 2n+li x ) =x i -^^ 2 F 1 (-^/V-f) • (22) 
When 7 2 = 1, the solution is of a different type and given by 



P2n + l(x)=x(- 1 rJ { I±^ 1 F 1 ( a n -y\. (2.1) 



n!/3 1 + l 

The proof is rather straightforward. It can be deduced from certain forward or backward shift 
operator formulas for the orthogonal polynomials that can be identified with the above expressions 
(see later). Alternatively, (|22p follows from the following contiguous relations for (terminating) 
hypergeometric series: 

s „ f—n.a \ , . „ f—n + l.a \ , „ /— n, a— 1 
(& + n) 2 iq( 6 + 1 ;^J -n(l-z) 2 F 1 { &+ 1 J = b 2 F X ( ' ft :: 

( —n,a \ (—n — \,a \ az f-n,a + l \ 
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Such contiguous relations are trivial to verify by comparing coefficients of z in left and right hand 
side. Similarly, (|23|) follows from: 

> „ / —n \ „ f—n + 1 \ „ /— n 

^(»-^(T^)=f'«G»- < 25 » 

Now it is a matter of identifying the above polynomials in order to find For this purpose, 

recall the definition of the Meixner polynomial M n (fc; /?, c) of degree n in fc, with parameters f3 and 

c [HEIEB]: 

M n (fe ;/ 9,c)= 2 F l (~ n, ~ k ;l--Y (26) 
These polynomials satisfy a discrete orthogonality relation: 

V^c fc M m (fe;/3, C )M n (fe;/3, C ) = ° ^ 5 mn (27) 

when /3 > and < c < 1. 

Whether the polynomials in Proposition [2] can be identified with Meixner polynomials, depends 
on 7. We should distinguish three cases. 
Case 1: (7! > 1. 

From (|22p it is clear that we can identify c with I/7 2 , and we have 



n! 7^ — 1 7^ 



p 2n+1 (x) = x (_ 7 )-« A /^±ikM n (^- - 1; /3 + 1, 4)- (28) 

n!p 7 Z — 1 7 Z 



The orthogonality relation (I27p leads to the following result: 

Proposition 3 For | 7 | > 1, the polynomials p n (x) satisfy a discrete orthogonality relation: 

/ 2 \ P 

w(x)p n (x)p m (x) = ( ? 7 - j <5 m „, (29) 

where 

Si = {±^7 2 - I fc E Z + }, (30) 
and where the weight function is given by 

w(x) = 1 for x = 0, 

^( x ) = I^k 7 - 2fc /or x = i^T^ly^ (fc = 1,2,...)- (31) 



So in this case, the spectrum of the position operator q is discrete and given by (|30p . 
Case 2: | 7 | = 1. 

For I7I = 1, we have already deduced in Proposition [2] that the polynomials are ±Fi series, and 
these can be identified with (generalized) Laguerre polynomials. 



P2n(x) = (- 7 )»y-^L(f-D(*»), p 2n+l {x) = (_ 7 )»y_^_ xL^(x 2 ). (32) 
The orthogonality relation of Laguerre polynomials leads to the following result: 
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Proposition 4 For \j\ = 1, the polynomials p n (x) satisfy a continuous orthogonality relation: 



L 



+00 

w{x)p n (x)p m {x)dx = r(P)5 mn , (33) 

00 



where 

w{x) = e~ x2 \x\ 213 - 1 . (34) 

So in the second case, the spectrum of the position operator q is continuous and given by M. 
Case 3: < | -y | < 1. 

Having found Meixner and Laguerre polynomials for the first and second case, one might expect 
to find Meixner-Pollaczek polynomials for the third case. However, this is not so. To see the proper 
form, one should first apply a transformation on the 2-F1 series in (122D. 



, -n,a \ , . _ „ / —n,b — a z , 
it\ ( b ;z)=(l-zr n 2 F 1 ( > b ;— ]. (35) 



After this, one can again identify the polynomials with Meixner polynomials (now with c = 7 



2V 



nl 1 — 7 Z 



P2n(x) = (-j) n \ M„(- . - /3; /3, 7 



2\ 



p 2n+1 (x) = x(- 7 ) n 1 / ^7^M n (- - /3; /3 + 1, 7 2 ). (36) 



n!/3 M «4-7 2 
The orthogonality relation (I27p now leads to the following result: 

Proposition 5 For < I7I < 1, the polynomials p n {x) satisfy a discrete orthogonality relation: 

^2 W(x)p n (x)p m (x) = — 2 ,p 5mn, (37) 

where 

S 3 = {±Vl-7V/? + ^ I fe G Z+}, (38) 

anc? where the weight function is given by 

W{x) = \^kt l2k forx = ±^/T^^pTk (A; = 0,1,2,...). (39) 

So in the third case, the spectrum of the position operator q is again discrete and given by (f38|) . 
Case 4: 7 = 0. 

In principle, there is a fourth case with 7 = 0, but this is somewhat trivial and we shall not 
return to it later. Indeed, for 7 = the matrix M q falls apart into irreducible (2 x 2)-blocks, 
because all S n = 0. The spectrum corresponds to the eigenvalues of these (2 x 2)-blocks, and it 
turns out that the corresponding polynomials are just discrete delta-functions: 

Proposition 6 For 7 = 0, the polynomials p n (x) satisfy the discrete orthogonality relation: 

w(x)p n (x)p m (x) = 5 mn , (40) 

where 



5 4 = {±y/f3 + k I k G Z+}, (41) 
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and where the weight function is constant: 

w(x) = - forx = ±y / P + k (k = 0,1,2,...)- (42) 

The polynomials are 

P2n(x) = <5 n ,fc for X = ±\f + k, 

P2n+i(x) = ±5 n)k for x = ±y/fi + k. (43) 

The above analysis determines the spectrum of the position operator q in the representation 
determined by f3 > 0, and for all possible values of the parameter 7 in (|16p . For each x belonging 
to the spectrum of q, the corresponding formal eigenvector is given by (12ip . Since essentially only 
7 2 plays a role, and since 7 = is a redundant case, we shall in the rest of the paper deal with 
7 > 0. 

Now the determination of the spectrum and eigenvectors of the momentum operator p is a 
formality. Due to the simple connection between (I16j) and fjlTf) . one can easily deduce that the 
spectrum of p is the same as that of q (in all four cases). Furthermore, the formal eigenvector of p 
for the eigenvalue y is given by 

00 

w(y) = J2(-i) n Pn(y)\P,n), (44) 

n=0 

where the p n 's are the same polynomials that appear in the analysis of q. 



4 On the shape of position and momentum wavefunctions 

The position (resp. momentum) wavefunctions of thes[(2|l) finite oscillator are the overlaps between 
the (normalized) g-eigenvectors (|2ip (resp. p-eigenvectors (|44p ) and the //-eigenvectors. Because 
of the close relation between (|21h and (j44]) it will be sufficient to study only the position wavefunc- 
tions. Obviously, these wavefunctions depend on the representation parameter /3 > 0, and on the 
parameter 7 which appears in the expression (I16p of q, so we will denote them by $^ (#): 

00 

v(x) = J2^\x) \M- (45) 

n=0 

Herein, x belongs to the spectrum of q, and $^f' 7 ^(x) is the polynomial p n (x) as in (f2Tj) but 
normalized. 

Let us begin with a familiar case, namely for 7 = 1, when the spectrum of q is R. From (|33|) 
and (I34p . it follows that the normalized versions of (I32p are given by 



^S(^) = (-W r(n+^ + l) l-l^ 1/2 ^ 2/2 ^(- 2 )- (46) 

These are just the paraboson wavefunctions, see e.g. [TUJ (A. 11)] (with /3 equal to the paraboson 
parameter a in [H]). In particular, when the paraboson parameter equals 1/2, one is just in the 
ordinary boson case; and indeed one has: 

d,a/2,D (x) = 1 *-* 2 / 2 HJx), (47) 
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with H n (x) the common Hermite polynomial. So for (/3, 7) = (1/2, 1), the sl(2|l) oscillator coincides 
with the canonical oscillator, whereas for 7) = 1) it coincides with the paraboson oscillator 
with paraboson parameter f3. 

Let us now consider the cases < 7 < 1 and 7 > 1. For < 7 < 1, the expressions of the 
(discrete) wavefunctions follow from (f36]) . using ([39]) and the normalization (|37"|) : 

*£ 7) (z) = (- X ) n W^%ir(l - 7 y/ 2 M n (A;;/?,7 2 ), 

*2#0»0 = (-l) n 7 B+ y ^^(! - W 2 *M n {k; P + 1, 7 2 ), (48) 

where s = ±\A = 7V/ 3 + fc > (fc = 0, 1,2, . . .)■ 
In a similar way, one obtains for 7 > 1: 

= (-l)^- n " fc ^7^^(7 2 - l)^M n (A; - l;/3 + 1, 1), (49) 

where x = ±vV- lVfe, (fc = 1, 2, . . .); 

for x = 0, the weight function is not simply the last expression in (|3ip . and therefore we have a 
separate expression for the wavefunction at x = 0: 



*£">(0) = (-1)^7— ^a/%7 2 " I)' 372 , $£?i(0) = 0. (oO) 

V • 

Note that we could also keep the first expression in (H9j) and multiply it by v/T-P <5fc,o; then it 

coincides with (|50p for = (or x = 0), and then we have a unified expression for all k = 0, 1, 2, 

Let us now consider the plots of these wavefunctions <E>^f ' 7 ^(x), for some values of the parameters 
j3 and 7 and for some n- values. As a first case, it is interesting to take /3 = 1/2, since we know 
that the case (/3,7) = (1/2,1) coincides with the canonical quantum oscillator. In Figure [T] we 
have plotted the 'ground state' wavefunction <£q 1//2,7 ^ (x) and the 'first excited state' wavefunction 
$S 1/2,7) (x), for some values of 7. For 7 = 1, this yields the common wavefunctions of the canonical 
oscillator, with support R. Then, we have plotted the wavefunctions for some values of 7 < 1, 
where the expression is used; and for some values of 7 > 1, where the expression (|4"9j) is used. 
In both of these cases, the support of the wavefunction is discrete (but infinite), so the plots consist 
of an infinite number of dots (of course, in the figures we can show only a finite number). Observe 
the similarity between the discrete plots of <&„ 1/ ' 2 ' 7 ^ (x) for 7 / 1 and the continuous wavefunction 



&n (%)■ When 7 tends to 1 (either from above or from below), the discrete plots of $n ( 
tend to the continuous plot of ^l 1 ^ 2 ' 1 ^ (x), provided the dots in the discrete plot are prope 
redistributed. This is necessary because in the continuous case the function satisfies 

j —00 



whereas in the discrete case one has 

2 



EK /2,7) K) =1 (7/1), 

xeS 
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with S the support given by (j30|) or (|38p . This observation also follows from the limits determined 
at the end of this section. 

In Figure[2]we have plotted the 'ground state' wavefunction <j?Q /3 ' 7 ' ) (x) and the 'first excited state' 
wavefunction $^' 7 ^(x), again for some values of 7, but now for another value of j3: /3 = 2. For 
7 = 1, this yields the wavefunctions of the paraboson oscillator, with support R. The wavefunctions 
for the other values of 7 yield discrete plots that tend to the paraboson wavefunctions when 7 tends 
to 1. 

Let us briefly return to the limits of the discrete wavefunctions when 7 — > 1. In the case 
< 7 < 1, consider the even wavefunction in (I48p . The essential limit comes from 

2 1 
limM„(fe;/3,7 2 ) = lim M n (-^ - /3; 0, 7 2 ) = y^-L^V), (51) 

7->l 7->l 1 - 7 Z (p)n 

which is a slightly modified form of a known limit [26, p. 243]. For the odd wavefunctions, the limit 
is similar, and also the factor coming from the weight function is easily computed under the limit, 
so one finds indeed that the expressions (I48h yield those of (I46h in the limit 7— >■ 1 (0 < 7 < 1): 

lim $(f' 7) (z) = 

The limit 7—7-1 for the case that 7 > 1, i.e. the wavefunctions ()49|) . is the same and the computation 
is similar to the one just described. 

In this section, we have paid attention only to the position wavefunctions. The momentum 
wavefunctions are completely analogous, and in fact it follows from (I44p that they are given by 

^Hy) = (-i) n ^(y), (52) 
where the last expression is that of the position wavefunction. 



5 Expressions for the sl(2|l) Fourier transform 

In canonical quantum mechanics, the momentum wavefunction (in L 2 (1R)) is given by the Fourier 
transform of the position wavefunction (and vice versa), with kernel K(x,y): 



K(x,y)$>(x)dx, K(x,y) 



1 



-ixy 



2vr 



This means that K{x, y) is also the overlap of the (formal) position eigenvector for the eigenvalue 
x with the momentum eigenvector for the eigenvalue y. 

So in the current case, the kernel for the corresponding s[(2|l) Fourier transform is given by 

00 00 

kV«\x, v ) = (v(x),w(y)) = X>^ 7) (*)«f' 7) (y) = E(-0"^ ,7) (^' 7) (z/)- (53) 



n=0 



n=0 



We need to compute this function in three cases, according to0<7< 1, 7 = 1 or 7 > 1. Let 
us start with the known case 7 = 1. The corresponding kernel has been computed in a number of 
papers [3T1I34]: 



K^\x,y) 



\xy 



13-1/2 



0F1 



x 2 y 2 



2/ 3 r(/3) 

It is easy to see that for j3 = 1/2 this expression reduces to 

1 1 
[cos(xy) — isin(xy)) 



x y pi 



x 2 y 2 



(54) 



-ixy 
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For the remaining cases, one can make use of the bilinear generating function for Meixner 
polynomials [3] (which can be deduced from [12], eq. (12), p. 85] or [391 proposition 3]): 



n=0 



^^M n (£;&,c)M n (£';M = (1 - z)- b -^'(l - £)W 2 F X 
n\ c 



b ' (z-c) 2 



(55) 



Using this, we obtained for the case < 7 < 1: 



K^\x,y) 



1 



27 



2 \\ + r 



1 + 7 2 



(P)Mi 



-k,-l 



4 7 



/3(i+7 2 : 



fcU! 

2-^1 



+ 1' 4 17 7 J 



(56) 



where x = ±\/l — j 2 \/ /3 + k, y 
In a similar way, one obtains for 7 > 1: 



± v / l- 1 2 ^/p + l, (k,l = 0,1,2,...) 



2 7 ^' £+ 77 2 -l^^/(/3) fc (/3) ^ 



-k,-l 

P 

where x = ±\/7 2 ~~ lVre, 



+ 7 2 

4 U 7 ; 



7 2 + l 



k\l\ 



. xy{l + j 2 ) 

1 TT^-s 2*1 



(i + <y fc , ) (1 + ^,0 

*! + 1,-Z + 1 1 



1 



4/3 7 2 " 1 V Z 3 + 1 
±V7 2 " 1 A (A:,/ = 0,1,2,...). 



t(7 "-) 2 
4 7 



(57) 



Note that one can explicitly compute the limits 7 — > 1 of the kernels (|56p and (|57p . and show 
that they yield the known paraboson kernel (|54j) . Such computations are based on limit relations 
of the following type: 



lim 2-^1 
7— >i 

7<1 



0-^1 



/3 ' 4 17 7 j 



lim 2-fi 
7— >1 

7<1 



/9 



>/3- 







■--(7--) 

' 4 17 7 j 



x 2 y 2 



6 The paraboson oscillator and osp(l|2) C b[(2| 1) 

Wigner |41| introduced the one-dimensional Wigner oscillator or paraboson oscillator [32J, leading 
to the field of Wigner quantization [23|. 133] . We shall recall some formulas for the paraboson 
oscillator (see [IS] or the appendix of |19j). In terms of the momentum and position operator p 
and q, the Hamiltonian of the paraboson oscillator is given by 



p 2 q 2 



«0 = * T + \. (58) 

Wigner dropped the canonical commutation relation [q,p] = i, but required instead the compati- 
bility between the Hamilton and the Heisenberg equations. These compatibility conditions are: 



p 2 q 2 p 2 q 2 



-ip. 



(59) 



So, one has to find self-adjoint operators p and q, acting in some Hilbert space, such that the 
compatibility conditions (|59p hold. The solutions to (|59p can be found by introducing two new 
operators b + and b~ (the paraboson creation and annihilation operators): 



V2 



{q Tip), 



(60) 
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or equivalently 

q = ±=(b + + b-), p=-L(b + -b-). 

It is then easily to see that Hq = , 6+}, and that the compatibility conditions (|59p are equivalent 
with 

[{6",6+},6+] = ±26+. (61) 

The relations (|61|) are nowadays the defining relations of a pair of paraboson operators 6+ |16j . 
Furthermore, it is known that the Lie super algebra generated by two odd elements 6+ subject to 
the restriction (|6ip is the Lie superalgebra osp(l|2) p3]. Keeping in mind the self-adjointness of q 
and p, i.e. (6+)^ = 6^, one is then faced with finding all star (or unitary) representations of the Lie 
superalgebra osp(l|2). These are known, and are characterized by a positive real number /3 and a 
vacuum vector |/3,0), such that 

6-^,0) = 0, {6~ 6+}|/3,0> = 2/3|/3,0). 

The representation space is the Hilbert space £ 2 (Z, + ) with orthonormal basis vectors |/3, n) 
(n € Z + ) and with the following actions: 

b + \f3,2n) = v / 2(n + /3)|/3,2n + l), 6~|/3,2n) = V2n|A2n- 1), 



(62) 



6+1/3, 2n + l) = v /2(n + l)|/3,2n + 2), 6"|/3,2n + l) = y^n + /3) |/3, 2n), 

from which follows 

{&-, 6+}|/3, n) =2(n + /3) |/3, n), (63) 

leading to the spectrum of Hq. From the action of Hq and from the explicit action of the commutator 
[q,p] on the basis vectors |/3,n), it is clear that for /3 = 1/2 the paraboson oscillator yields the 
canonical oscillator. The position wavefunctions for the paraboson oscillator can then be determined 
by constructing the formal eigenvectors of q (see [19]) or by different techniques |31j . 

To see that the Lie superalgebra generated by the paraboson operators is indeed osp(l|2), let 

H= 1 -{b-,b + }, £+ = i{6+,6+}, E~ = -\{b~, 6-}; B+ = ^=6+, B~ = ±=b~ . (64) 

Then, using only (f6T|) one finds back the standard osp(l|2) commutation relations [131 p. 260]: 

[H, E±] = ±E ± , [E+,E~] = 2H, [H, B±] = ±^B ± , 

[ J B ± , J B+] = - J B ± , {B+,B~} = ~H, {B±,B ± } = ±~E ± . (65) 

From these relations, the embedding of osp(l|2) into sl(2|l) is also clear. Starting from ([4])-([6]), 
keeping H, E + and E~ and putting 

B + = l -(F + + G + ), B- = \{F- - CT), (66) 

one finds again the relations (|65p (which is why we have used the same names for the corresponding 
generators) . 

It is now easy to verify that the sl(2|l) irreducible representation ILj of section 2 decomposes, 
under the embedding sl(2|l) D osp(l|2), as a single irreducible representation Tp of osp(l|2). Note 
that in this context the paraboson position operator q is given by 

q = ^(6+ + b~) = 2(B + + B~) = F + + G + + F~ - G~ . (67) 
v2 

Comparing with (|16|) reveals why the case 7 = 1 corresponds to the paraboson oscillator, and why 
the st(2|l) oscillator wavefunctions $^'^(a;) coincide with the paraboson oscillator wavefunctions. 
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7 Remarks and discussion 



The oscillator models described by the discrete series representations Ilg of the Lie superalgebra 
s[(2|l) offer an interesting extension of classical oscillator models. They are characterized by two 
parameters: /3 > is a representation parameter, and 7(7^ 0) is an extra parameter appearing in 
the expression for the position operator q. 

Let us consider here some other quantities in these models that may play a role in physics. First 
of all, using the expression of q and p and the actions (|10p . one finds in the representation 11^: 

[q,p] |/3, 2n) = 2i (/3 + (1 - 7 2 )n) |/3, 2n) (n = 0, 1, 2, . . .), 

|/3, 2n-l) = 2i ((1-/3)- (l- 7 2 )n) |/3, 2n-l) (n = 1, 2, . . .). (68) 

So, just as for the paraboson oscillator, the action of the commutator [q,p] is still diagonal. Note 
that for 7 2 = 1, one finds indeed 

[q,p] |/3, 2n> = 2i/3|/3,2n> (n = 0, 1, 2, . . .), 

[q,p] |/3, 2n - 1) = 2i(l - /3)|/3, 2n - 1) (n = 1, 2, . . .), 

which is a known expression for the paraboson case. And clearly, for (3 = 1/2, this becomes 
[q,p] = i, the canonical situation. 

Another operator that is worth considering in 11^ is \ + \, since this stands for the Hamiltonian 
in the paraboson (and thus also in the canonical) case. It is easy to verify that 

p 2 q 2 /7 2 + 1 \ 

{— + — )|/3,n) = I — - — n + /3J|/3,n) forn even, 
V 2 q 2 fl 2 + l 7 2 -l\ 

( y + n) = [ 1 ^ n + P + -^-J I A n ) fOT « odd - ( 69 ) 

Again, this operator is diagonal in the general case, with a very simple action. Clearly, for 7 2 = 1 
this reduces to the paraboson oscillator Hamiltonian (or the Wigner oscillator). 

To summarize, we have developed a new model for the quantum oscillator based upon the 
Lie superalgebra st(2|l) and its discrete series representations n^. These are infinite-dimensional 
unitary representations labeled by a positive number /3, and the action of the sl(2|l) basis elements 
in this representation has been determined explicitly. The Hamiltonian H, the position q and 
the momentum p of the model are three self-adjoint elements of sl(2|l) satisfying the Hamilton- 
Lie equations (JTJ) . In particular, this requirement leaves a degree of freedom in the choice of the 
position operator, giving rise to an arbitrary parameter 7 in the expression for q. The spectrum 
of H coincides with that of the canonical quantum oscillator. The spectrum of q depends on 7, 
and can be infinite discrete (I7I 7^ 1) or continuous (| 7 | = 1). The position wavefunctions ^n'^\x) 
have been determined explicitly. For 7 = 1 they coincide with paraboson wavefunctions, and are 
given in terms of Laguerre polynomials. In particular, the wavefunctions $^ 1/ ' 2 ' 1 ' ) (x) are those of 
the canonical quantum oscillator in terms of Hermite polynomials. For | 7 | 7^ 1 they are given in 
terms of Meixner polynomials, and satisfy a discrete orthogonality relation. Plots of the discrete 
wavefunctions reveal properties that are very similar to those of the canonical oscillator (when 
(3 = 1/2) or to the paraboson oscillator (when /3 > 1/2). From the closely related momentum 
wavefunctions, the sl(2|l) Fourier transform has been constructed in explicit form. The embedding 
of the Lie superalgebra osp(l|2) into s[(2|l) offers an algebraic explanation of the appearance of 
paraboson wavefunctions in the sl(2|l) model, since osp(l|2) is the superalgebra underlying the 
paraboson oscillator algebra. 
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The wavefunctions appearing here are Meixner polynomials of the type M n (x 2 ) for even wave- 
functions and of the type xM n (x 2 ) for odd wavefunctions. This kind of structure is reminiscent of 
some so-called "—1 polynomials". The —1 polynomials are usually considered [40J to be the q = — 1 
limit of basic hypergeometric polynomials appearing in the g-Askey scheme [26J (as long as this 
limit makes sense). As far as we know, they have not been explored systematically [3D]. The first 
example appeared in work of Bannai and Ito [9], where a q = — 1 limit of g-Racah polynomials 
plays a role. More recently, the q = — 1 limit of little g-Jacobi polynomials was considered in |40| . 
and the q = — 1 limit of dual g-Hahn polynomials in }38] , For these examples, expressing the — 1 
polynomials in terms of their classical counterparts P n (x), one sees indeed that they are essentially 
given in terms of P n (x 2 ) and xP n {x 2 ). 

Note also that the underlying algebra sl(2|l) has a natural coproduct, so it should be feasible 
to construct tensor products of the representations given in this paper, and determine the corre- 
sponding Clebsch-Gordan coefficients. Note that the representations IL3 of sl(2|l) coincide (see 
section 6) with the representations Tp of osp(l|2), as far as the representation space and the action 
of osp(l|2) is concerned. These paraboson oscillator representations are essentially equal to the 
class of representations of sL_i(2), considered in [15]. The Clebsch-Gordan problem for s[_i(2) has 
been solved in |15j . and the Clebsch-Gordan coefficients are determined in terms of dual —1 Hahn 
polynomials. It would be interesting to see whether the Clebsch-Gordan problem for sl(2|l) has the 
same solution. If that is the case, this might lead to interesting relations between dual —1 Hahn 
polynomials and the Meixner polynomials appearing in this paper. 
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Figure 1: Plots of the wavefunctions <J>^ (x) in the representation with /3 = 1/2, for n = (left 
column), and n = 1 (right column), for 7 = 0.4,0.75,1.0,1.2,1.5. For 7 = 1 the wavefunction is 
continuous, for the other 7-values it is discrete. 
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Figure 2: Plots of the wavefunctions $n \%) in the representation with /3 = 2, for n = (left 
column), and n = 1 (right column), for 7 = 0.4,0.75,1.0,1.2,1.5. For 7 = 1 the wavefunction is 
continuous, for the other 7- values it is discrete. 



19 



